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An  earlier  paper  reported  a  theoretical  study  of  the  ettects  ot  torouial  forces  in  current  loops  embedded 
in  an  ambient  plasma  such  as  the  corona.  Ihese  forces  occur  m  current-carrv  mg  plasmas  with  curvature 

The  model  loop  is  characterized  by  twisted  magnetic  field  lines  1  he  magnetic  held  and  current  densitx  have 
toroidal  and  potoidul  components.  In  this  earlier  worV.va  current  loop  was  assumed  to  evolve  while  maintain¬ 
ing  a  half-torus  configuration.  The  important  geometrical  constraint  that  tootpomts  are  general  Iv  immobile  on 
the  relevant  time  scales  was  not  taken  imo  account  In  the  present  paper,  the  dv  nautical  properties  of  current 
loops  similar  to  those  of  the  earlier  work  are  examined  with  the  requirement  that  the  tootpomts  remain^sia- 
tionarv  in  the  photosphere  Other  essential  ingredients  ot  the  model  are  unchanged  to -t* --found  -that  .ibe 
behavior  of  the  present,  more  realistic  nit  Kiel  is  qualitativelv  similar  to  that  of  the  earlier  model  fhis  is  due 
to  the  fact  that  the  toroidal  forces  have  only  a  mild,  logarithmic  dependence  on  the  local  aspect  ratio  In  par 
Ocular,  the  velocities  which  the  apex  can  achieve  under  the  action  of  toroidal  forces  and  the  magnetic.  energx 
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released  in  the  form  of  drag  heating  are  similar  in  iheir  respective  magnitudes.  In  the  present  calculation,  the 
ratio  r}  of  the  height  of  the  apex  and  the  tooipotm  separatum  is  a  physical  parameter  It  is  tound  that  a  tow 
lying  loop  (small  q)  ean  evolve  to  a  "taller'*  (larger  ij>  loop  in  a  quasi-equilibrium  tushion  it  the  current  is 
increased  slowly  while  holding  the  footpoint  separation  fixed  Tune  evolution  the  loops  and  the  magnetic 
energy  converted  via  drag  heating  are  presented  Results  are  also  presented  tor  loops  with  relatively  strong 
current.  The  results  are  discussed  in  the  context  ot  the  solar  environment.  Itt-pttrtKHditr4;  it  appears  that  the 
toroidal  forces  can  lead  to  a  wide  range  ot  plasma  motion.  The  a»vomp;»nving  drag  heating  can  take  place 
with  a  eorres fn  Hid i.iglv  ...uc  :a  :  •;  of  rates 
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DYNAMICAL  PROPERTIES  OF  SOLAR  CURRENT  LOOPS 
WITH  LINE-TIED  FOOTPRINTS:  EFFECTS  OF 
TOROIDAL  FORCES 

INTRODUCTION 


In  an  earlier  paper  (Chen,  1987),  the  dynamical  behavior  of  a  current 
loop  under  the  action  of  curvature  ("toroidal")  forces  was  studied.  The 
model  loop  is  immersed  in  an  ambient  plasma  such  as  the  corona  and  is 
characterized  by  twisted  magnetic  field  lines.  The  magnetic  field  and 
current  density  have  toroidal  and  poloidal  components.  The  model  assumes 
that  the  current  conservation  is  established  in  or  below  th<_  photosphere. 
The  calculation  showed  that  if  the  apex  position  of  an  equilibiium  loop  is 
perturbed,  the  loop  can  expand.  The  condition  for  instability  to  such 


expansion  was  given  in  terms  of  a 
magnetic  flux  enclosed  by  tile  loop 


parameter  e  =  where 

above  the  photosphere  and 


<f> 

P 

4> 

T 


is  the 
i  s  t  he 


total  flux  of  the  entire  current  structure  including  the  subpho t ospher i c 


current.  Although  the  subphotospheric  flux  or  current  structure  is  not 


measurable  in  reality.  it  is  a  physically  meaningful  quantity.  Tire 
analysis  used  an  integral  form  of  f  c x  B  Vp  to  describe  the  centei 
of -mass  motion  of  the  apex.  The  drag  between  expanding  Loops  and  ambient 
gas  was  included  using  a  simple  model.  This  allows  release  of  magnetic 
energy  via  drag  heating.  A  less  developed  predecessor  of  this  model  was 
discussed  by  Xue  and  Chen  (1°80)  in  the  context  of  a  magnetic  energy 
release  mechanism. 


The  main  results  arc  the  following.  (1)  There  exists  a  critical  value 

E  which  depends  only  on  physical  parameters  above  the  photosphere  .such 

that  an  equilibrium  loop  with  r  <  e  (r  >  ef.  )  is  unstable  (stable).  (?) 

An  unstable  loop  can  expand  and  sometimes  teach  a  "second"  equilibrium 

following  a  period  of  damped  small  oscillation.  A  loop  initially  in 

equilibrium  typically  expands  at  subsonic  velocities.  This  corresponds  to 

slow  drag  heating  of  the  ambient  gas.  (3)  For  a  loop  which  is  not 

initially  in  equilibrium  (perhaps  due  to  loss  of  equilibrium)  with  moderate 

to  large  current  and  magnetic  fields,  the  apex  can  bo  driven  super  son  i  ra 1 1 y 

s 

through  the  ambient  gas.  For  a  loop  with  ~30r;  initially  at  a  height  of  It' 

km,  the  aoex  can  attain  velocities  ol  loughlv  1300  km  s  1  foi  tens  of 

.32 

minutes,  dissipating  tip  to  10  erg  via  diag  heating  (also  shock  heating  in 
the  supersonic  case).  During  the  expansion  ol  loops,  the  l.orentz  force  ran 


do  work  so  that  the  magnetic  energy  ran  be  dissipated.  Thus,  a  vide  range 
of  apex  velocities  and  magnetic  energy  release  rates  are  possible  within 
the  context  of  the  model. 
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Curvature  forces  are  a  particular  iorm  of  Lorentz  and  pressuie  forces 
which  occur  in  curved  segments  ot  current-carrying  plasmas  and  act  along 
the  local  radius  of  curvature.  (A  similar  force  also  occurs  in  metallic 
current  hoops.)  In  the  laboratory,  these  forces  are  wel 1  -  under s t ood . 
However,  laboratory  plasmas  are  typically  surrounded  by  vacuum  which  in 
turn  is  enclosed  in  rigid  metallic  containers.  In  addition,  magnetic 
fields  are  applied  by  external  coils  to  balance  the  curvature  forces 
(sometimes  referred  to  as  the  "hoop  stress").  In  the  solar  and 

astrophysical  environments,  magnetic  structures  are  usually  embedded  in 
plasmas  and  are  not  surrounded  bv  metallic  containers.  The  curvature 
forces  in  such  environments  have  not  been  fully  investigated. 

The  model  described  above  was  simplified  in  order  to  elucidate  the 

basic  physics  underlying  curvature  forces  and  their  effects  in  solar  and 

astrophysical  environments.  One  simplification  is  that  the  initial 

geometry  of  the  loop  is  a  half  torus  <>l  uniform  aspect  ratio  R/a,  where  R 

is  the  major  radius  and  a  is  the  minor  radius  and  that  the  loop  expands  (if 

unstable)  maintaining  the  half -torus  geometry.  The  uniform  aspect  ratio 

approximation  is  not  expected  to  result  in  a  serious  error  since  the 

curvature  forces  depend  on  the  local  aspect  tat.io  R/a  as  ln(8R/a),  a 

relatively  mild  dependence.  However.  in  a  realistic  current  loop,  the 

footpoints  are  expected  to  Ire  essentially  immobile  compared  with  the  apex 

ort  the  relevant  time  scale  of  tens  ot  minutes  (the  so-called  photospheric 

"line  tying")  and  the  loop  does  not  remain  a  half  torus.  Tn  the  ideal  MHD 

approximation,  the  foot  point  magnetic  field  is  tied  to  the  "infinitely" 

conducting  photospheric  (and  subpho tospher i  <: )  plasma  which  is  much  denser 

than  the  coronal  plasma.  In  this  paper .  we  achieve  stationary  footpoints 

without  invoking  the  ideal  MHD  line  tying  because  of  the  small  fractional 

-  1  -4 

ionization  in  the  photosphere  (perhaps  10  "  to  10  ).  However,  the  ionized 

component  of  the  photospheric  plasma  is  coupled  to  the  magnetic  field  lines 


that  enter  the  photosphere.  The  ionized  particles  in  turn  are  coupled  to 

the  neutrals  via  collisions  on  the  relevant  time  scale.  If.  tire  magnetic 

flux  structure  extends  deeper  down.  the  fractional  ionization  is  expected 

to  increase.  Thus,  any  motion  of  the  footpoints  would  be  resisted  by  drag 

? 

forces  which  scale  as  n  (eq.  (  ID | )  where  n  is  the  photospheric  density 
and  Vj  is  the  footpoint  velocity.  Since  n  is  much  greater  than  the 
coronal  density  n  ,  the  footpoint  motion  due  to  similar  forces  m"s*  he 

5/2 

slower  by  (n^/n  )  in  comparison  with  possible  motions  of  the  apex.  if 

wc  estimate  n  ~  10^  cm  ^  and  n  ~  10^  cm  then  (n  /n  )  ^  ^  ~  If)  ^ . 
c  p  c  p 


In  the  present  paper,  we  will  use  a  simple  model  geometry  that 
explicitly  takes  into  account  the  immobile  footpoints.  The  basic  physics 
is  the  same  as  that  of  the  earlier  paper  (Chen,  1987).  Because  the  more 
realistic  geometry  requires  extensive  modification  of  equations,  it  seems 
desirable  to  document  the  calculation  in  detail.  We  therefore  repeat  the 
analysis  with  immobile  footpoints.  However,  the  interested  reader  is 
referred  10  the  earlier  paper  for  a  more  general  discussion  of  curvature 
forces  in  the  solar  environment  and  a  more  comprehensive  reference  list. 

We  will  start  with  a  model  current  loop  which  is  initially  in 
equilibrium  and  calculate  its  time-dependent  behavior  in  response  to 
perturbations  of  the  apex  height  (Sec.  II).  The  theoretical  framework  will 
be  first  presented,  followed  by  a  numerical  calculation  of  the  long-time 
evolution  of  loops  including  the  drag  force  due  to  the  ambient  gas  (Sec. 
III).  Curvature  forces  are  not  limited  to  equilibrium  loops.  As  an 
example,  we  discuss  the  behavior  of  a  loop  carrying  a  relatively  large 
current,  which  may  not  tie  in  equilibrium  initially.  Although  no  attempt  to 
model  specific  observation  will  be  made,  we  will  discuss  the  potential 
relevance  of  the  results  to  plasma  activities  in  t he  corona  (Sec.  IV).  It 
will  be  shown  that  a  cut  rent  loop  acting  under  the  influence  oi  cur  vature 
forces  can  exhibit  a  range  of  behavior  (expansion  velocity,  energy  release, 
etc.)  compatible  with  certain  motion  related  effects  in  the  solar 
envi ronmen  t . 


II.  DYNAMICS  OF  A  MODEL  CURRENT  LOOP 


In  this  section,  we  consider  the  evolution  of  an  isolated  loop  with 
twisted  magnetic  field  lines  which  is  initially  in  equilibrium.  Figure  1 
shows  schematically  a  model  loop  which  has  a  toroidal  magnetic  field  and 
poloidal  magnetic  field  with  a  poloidal  and  toroidal  current  density 
components  J  and  J  ^ ,  respectively.  The  apex  of  the  loop  is  at  a  height  of 
Z  from  the  photosphere,  R  is  the  local  radius  of  curvature  (major  radius) 
and  a  is  the  minor  radius.  The  aspect  ratio  R/a  is  taken  to  be  roughly  5 
to  10.  The  loop  is  embedded  in  a  field-free  plasma  of  pressure  p  .  In 

3 

order  to  satisfy  current  conservation,  we  allow  the  current  to  close  in  or 
below  the  photosphere  where  the  plasma  is  simply  assumed  to  be  much  denser 
than  the  coronal  plasmas.  The  footpoints  are  separated  by  a  distance  2so 
and  are  assumed  to  be  essentially  immobile  because  of  the  dense 
subphotospheric  plasmas.  No  particular  current  distribution  will  be 
specified  below  the  photosphere. 

The  ambient  plasma  above  the  photosphere  is  assumed  to  decrease 
exponentially  with  a  scale  height  H.  In  the  corona,  H  can  be  given  by 


2kT 


where  k  is  the  Boltzmnn  constant.  T  is  the  ambient  plasma  temperature,  m. 

a  4  1 

is  the.'  ion  mass  and  g  is  the  gravitational  acceleration  which  is  2.7x10  cm 

-  2  5 

sec  at  the  surface.  At  the  base  of  the  corona,  H  is  roughly  10  km. 

A.  Curvature  forces 

Curvature  forces  can  occur  in  any  curved  segments  of  current-carrying 
plasmas  and  depend  on  the  local  major  radius  R  and  local  minor  radius  a. 
In  a  solar  loop  configuration,  we  expect  the  apex  region,  say,  one  third  of 
the  loop  about  the  apex,  to  exhibit  t  Ik  greatest  degree  of  motion  because 
it  is  farthest  from  the  neatly  immobile  footpoints.  The  motion  of  the  apex 
region,  schematically  shown  in  Figute  1  as  the  region  between  the  two 
dashed  lines,  is  determined  by  the  Loientx  foice,  which  consists  of  both 
the  magnetic  tension  and  pressnic  fotces,  and  Vp .  The  local  foice  density 
f  acting  on  a  plasma  element  is  gi  < n  by 

f  f-  J  x  B  Vp,  (la) 
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with  the  usual  equation  relating  the  cut  rent  and  magnetic  field 

J  =  -^  V  x  B .  (lb) 

—  4  n.  —  — 

In  this  paper,  the  displacement  current  is  neglected.  Each  plasma  element 
moves  according  to  £.  However,  the  motion  of  the  center  of  mass  of  each 
segment  of  the  loop  is  determined  by  the  force  integrated  over  the  given 
segment.  In  considering  the  motion  of  the  loop  apex,  we  will  refer  to  this 
center-of-mass  motion.  The  center  of  mass  of  a  slice  of  t  lie  loop  lies  on 

the  dash-dot  line  along  the  loop.  Tn  order  to  integrate  f  over  segments  of 

the  loop,  we  make  a  geometrical  simplification  that  t he  loop  above  the 
photosphere  is  a  part  of  r  toroidal  plasma  which  is  intersected  by  the 
photosphere  with  a  footpoint  separation  of  2s  .  The  current  in  oi  below 
the  photosphere,  which  conserves  current,  will  not  be  specified  in  detail. 
For  such  configurations,  we  see  from  Biot-Savai  t  law  that  the  contribution 
to  the  magnetic  field  neat  the  apex  due  to  the  cut  tent  segment  below  the 
photosphere  scales  as  (a/R)0'n.  vhete  0  is  the  angle  between  the  vertical 
and  the  line  cwinP'’  t  i  ng  a  footpoint  to  the  axis  of  the  torus  (see  Fig.  1 
and  the  definition  of  8  lol lowing  eq.  [  1 1  J ) .  In  view  of  this 

simplification,  we  will  generally  consider  loops  which  arc  large)  than 

half-torus  fot  a  given  secant  ?$  .  Then ,  0/ n  is  less  than  If.  Because  ve 

consider  R/a  in  the  range  of  5  to  in,  a.R  <<  1.  As  the  apex  height 
increases,  the  loop  becomes  moi e  toroidal  with  decreasing  0/ n  <<  1 . 
reducing  the  correction  term.  The  tieatment  is  expected  to  he  most 
accurate  for  curvature  forces  near  the  apex  which  terrains  nearly  semi 
toroidal.  This  approximation  is  similat  to  that  used  by  Anzer  (1^78)  and 
Van  Tend  (  1979).  This  geometry  is  an  over  -  simp] i f i cat i on  hut  it  has  the 
advantage  that  we  can  analytically  calculate  the  forces  acting  on  segments 
of  the  loop  near  the  apex  and  include  the  impoitant  feat  tit e  that  the 
footpoints  are  immobile.  Foi  an  imptoved  geometiv,  see.  f  <>t  example,  Auxei 
and  Poland  (1979).  Note  that  the  intpgi  a  1  lot'::  of  t  tale1  into  ,e  <  mint  the 
curved  geometry  in  a  naiuial  way. 

In  order  to  determine  the  motion  of  tin  apex  (  i  .  <;  .  ,  its  rentei  of 
mass),  we  integrate  f  over  a  section  of  the  touts  at  the  apex.  in  doing 
so,  we  must  impose  certain  conditions:  the-  loop  i-:  mim  minded  by  a  plasma 
of  pressure  p  (no  metallic  containet  suitouuding  the  loop)  and  the 
magnetic  field  vanishes  at  infinity.  A  technique  foi  integrating  ovet 
toroidal  plasmas  was  developed  previously  (Shafranov,  1H66),  originally  foi 


application  to  laboratory  plasmas,  and  the  analysis  will  not  be  reproduced 
here.  Using  the  same  analysis  hut  with  the  above  conditions  (as  opposed  to 
the  laboratory  conditions),  integration  of  1  over  a  section  of  the  torus 
yields  the  equation  of  motion  for  the  center  of  mass  of  the  apex 


where  F0  =  Md^Z/dt^  with  M 
K 

along  the  major  radius  per  unit 
the  minor  radius  a  also  changes. 


1  g  i  _1  t 

2  P  2  g2  ’ 

P 

2- 

rta  nm .  .  Here  F_  is  the  integrated 
1  K 

length  of  the  loop.  As  the  apex 

This  evolution  can  be  described  by 


(2) 


force 
moves . 


2  3- 
nc  a  nm. 


(3) 


Note  that,  for  the  simplified  geometry,  the  major  radius  R  is  related  to 
the  height  of  the  apex  Z  by 


R 


l. 

s 

o 


(4) 


This  equation  expresses  the  constraint  that  the  footpoints  are  immobile. 
Here  I  is  the  t  t a  1  toroidal  current  defined  by 


a 

2nj  dr  i J  . 
0 


i  ’ 

B  /Bn 
P 

where  p  is  the  average  interna]  pressure  of  the  loop,  p  is  the  ambient 
pressure  and  B  -  B  }(a)  is  the  poloidal  magnetic  field  at  the  outer  edge  of 
the  loop  (r  -  a).  The  quantity  £.  is  the  internal  inductance, 

characterizing  the  minor  radial  current  distribution,  and  ranges  from  0 
for  a  surface  distribution  to  1/2  for  a  uniform  current  distribution.  Note 
that  the  curvature  effects  are  relatively  insensitive  to  the  assumption  of 
uniform  R/a  because  of  the  logarithmic  dependence.  In  equation  (2),  mass 
flow  along  the  loop,  orthogonal  to  major  radial  expansion,  is  also 


The  quantity  £  is  defined  by 
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neglected  because  the  curvatuie  forces  occur  with  or  without  such  flow. 
Moreover,  mass  flow  is  important  only  if  the  flow  velocity  is  comparable  to 
the  Alfven  speed  in  the  loop. 

By  choosing  to  use  the  integrated  form  of  the  force  equations,  we 
sacrifice  the  description  of  each  plasma  element.  However,  this  is  not  a 
serious  disadvantage.  In  fact,  if  we  were  to  obtain  a  local  solution  for 
each  plasma  element,  we  would  integrate  the  result  to  determine  quantities 
such  as  the  center- of -mass  motion  in  order  to  compare  with  observation. 
The  integrated  result  would  then  take  on  the  form  of  equations  (2)  and  ( i ) . 
The  detailed  structure  enters  through  quantities  such  as  p  and  .  This 
technique  allows  one  to  obtain  the  integrated  form  without  first  finding  a 
detailed  point-wise  solution.  The  fact  that  MHD  is  amenable  to  global 
integrated  represen ta t i ons  has  proved  useful  for  various  problems  (e.g., 
energy  principles,  virial  tieoiem.  etc). 

In  Se~.  Ill,  we  will  integrate  these  equations  numeiieally.  In  this 
section,  we  will  first  use  a  simplified  version,  applicable  to  neai 
equilibrium  cases,  in  order  to  study  the  linear  behavior  and  gain  insight 
foi  the  later  results.  Note  that  if  the  loop  is  neatly  in  equ  i  1  i  hi  i  uni . 

H  n 

fotce  balance  along  the  minoi  tad  ins.  da'  '  d  t '  n,  gives 


Substituting  this  expression 


O 


1  ) 

c  I  B'  /B'  . 
t  p 

into  equation  (2),  we  obtain 
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which  describes  the  force  pet  unit  length  of  the  loop  acting  on  the  center 
of  mass  if  the  loop  is  near  equilibrium.  As  the  initial  configuration,  we 
will  adopt  a  model  loop  of  the  type  discussed  in  Xue  and  Then  (19B1.  For 
this  class  of  equilibrium  loops,  the  cutvatuie  forces  ate  explicitly 
balanced.  Here,  we  give  a  luief  summaty  of  equilihiium  proper  t  it's .  In 
equilibrium,  the  force  density  f  acting  on  each  element  of  the  loop  is 
zero.  Therefore,  we  have  t)  e:-:a<  t  1  y -  This  equat  ion  then  gives 


a 


p 
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Since  £. /2  is  generally  the  smallest  tetm,  ve  will  adopt,  lot  convenience', 
a  sui  face  current  model  and  set  -  0  hencefoi  th.  For  R/a  of  the  order  of 
10,  we  see  that  (3  <  0  in  equilibrium.  This  particular  equilibrium 
condition  arises  from  the  absence  of  ambient  magnetic  field.  This  point 
and  some  modifications  due  to  the  presence  of  ambient  fields  will  be 
discussed  briefly  in  Sec.  IV.  However,  the  essential  physics  of  cuivatute 
forces  is  not  limited  to  this  class  of  equilibria  or  by  this  condition.  It 
is  convenient  to  define  the  total  poioidal  current  by 


I 

P 


where  the  integration  is  over  the  minor  radius.  Fot  the  surface  current 
model,  we  have 


B 

t 


4  n 


1 

P 


with  B  -  t)  and  p  p  inside  the  loop.  tun  side  the  loop,  ve  have 


with  B.  -  it  and  p  u  and  i  i1'  measuted  along  the  minoi  radius  front  tiio 
t  1  f  a 

rente)  of  the  circular  rioss  section.  In  the  above  expressions,  the 
cot  recti  on  trims  of  the  ordei  (a'R)9/rt  due  to  the  geometric,-,! 
.simplification  are  neglected. 

The  above  expressions  ate  appropriate  for  current  carrying  plasmas 
embedded  in  an  ambient  plasma  with  no  metallic  containers.  As  noted 
before,  this  is  an  important  difference  from  such  laboratory  systems  as 
t okamaks .  Incidentally,  this  glass  of  equilibria  properly  satisfies  the 
requirements  of  the  virial  theorem  (Shafianov  19(i6). 


B.  Dynamical  Instability 

In  this  section,  we  investigate  the  stability  properties  of  the 
equilibrium  loop  with  tespect  to  pei  r  in  ha  t  ions  of  the  apex  height.  In  Sec. 
ITT.  we  will  consider  more  general  cases  with  Target  velocities.  In  the 
present  paper,  we  will  assume  lot  simplicity  that  R/a  is  uniform.  Ve  note 
that  the  curvature  forces  depend  on  the  local  aspect  ratio  as  ln(8R/a). 
which  is  a  rather  mild  dependence,  so  that  the  essential  physics  of 


H 


curvature  fotces  should  not  depend  sensitively  on  this  •approximation.  As 
the  apex  is  displaced  from  its  initial  equilibrium,  the  forces  experienced 
by  the  apex  region  be  given  by  linearizing  equation  (7)  which  is  valid  for 
small  deviation  away  from  equilibrium; 


d2(&Z) 

0 

d  t " 


1 1  CL  dR  Z  da  dgpl  6Z 

?ulI  IR  dZ  "  a  dZ  +  /jdZ  j  Z 

C  MK 


(8) 


2  . 

Mere  H  --  na  p  is  the  mass  pei  unit  length  of  the  loop,  6a  is  the  change  in 
the  minot  radius  and  p  is  the  average  mass  density  inside  the  loop.  The 


quail  t  i  t  y 


is  obtained  f l om  eq.  (5): 


(  If) 


whole  H  is  the  gravitational  sr.  ale  height.  Because  we  assume  that  t  hoi  e  is 
no  mass  flow  along  the  loop,  the  total  mass  M.f.  oi  the  loop  above  the 
photosphete  is  constant  in  time  where 


M  -  2  JT0RM  , 


wliei  e  9  =  sin  (s  /K)  and  ?rt9R  is  the  length  of  the  loop  above  the 
photosphere.  In  this  papei  ,  we  will  mainly  use  7.  >  s  .  Because  ve  use 
the  integiated  solution  rather  than  the  local  ( d i ( f ei en I i a  1  )  solution,  vr 
will  use  a  numhei  of  global  conditions  to  consttain  the  dynamical  hell. iv  ini  . 
Ue  assume  that,  oil  the  time  scale  of  the  loop  evolution,  possible  changt' 
in  the  fluxes  are  small.  Foi'  the  toroidal  flux,  we  assume 


9 


cons  tan  t 


(12) 


Here,  is  the  total  poloidal  flux  and  LT  is  the  total  self-inductance  of 
the  current  distribution  including  the  submerged  part  (Figure  1).  Note 
that  current  conservation  requires  only  that  there  be  some  current.  Given 
a  current,  the  total  flux  and  inductance  can  be  unambiguously  defined 
(albeit  not  necessarily  measurable)  without  specifying  details  of  the 
underlying  current  structure.  The  submerged  structure  is  included  in  the 


calculation  via  the  inductance.  Ue  can  define  the  inductance 
with  the  poloidal  flux  above  the  photosphere  by 


L  associated 
P 


L  = 
P 


$ 

_P 


where  the  total  poloidal  flux  is  $T  4> 


$  .  Then,  we  define  at  time  t 

s 


(14) 


This  quantity  c  is  a  rough  measute  of  the  relative  "size"  of  the  loop  above 
the  photosphere  and  the  entile  current  structure.  If  the  current  is  closed 
near  the  photosphere,  then  e  =  1.  If  the  submerged  current  structure  is 

much  larger  than  the  loop  above,  then  t  <<  1.  Note  that  the  initial  flux 
J  ,  thus  the  initial  inductance  of  the  submerged  current,  is  not 
calculated.  It  is  used  as  a  parameter  to  characterize  the  submerged 

current  structure.  Note  also  that  we  do  not  require  $  and  <1  to  he 

p  s  ' 

separately  constant  since  flux  emeigence  appeals  to  be  common. 

We  have  described  the  essential  ingredients  of  the  model.  We  will  now 
attempt  to  calculate  more  specific  properties.  For  the  dynamics  oi  the 
loop  interior,  we  assume  that  t  lie  adiabatic  expansion  law  is  valid: 


—  y  ~ 

p  V  =  constant 


where  y  is  the  adiabatic  index  and  where  V 


2  2 

2tt'a^0R  is  the  volume  of  the 
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loop.  Then,  we  have 


6p  = 


YP 
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From  equations  (4)  and  (11),  we  obtain 


and 
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Next,  from  toroidal  flux  conservation,  equation  (12),  we  obtain 

SB . 
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Sa 

a 


(16) 


From  the  definition  of  B  ,  we  find 

P 


SB 

_E.  _ 
B 

P 


(17) 


In  calculating  61^,  we  assume  that  the  submerged  current  structure, 
whatever  it  is,  remains  unchanged  on  the  relevant  time  scale  in  the  much 
denser  plasma.  This  assumption  is  consistent  with  equation  (4),  the 
immobile  footpoitits.  Then,  the  changes  in  L  are  primarily  due  to  changes 
in  the  loop  above  the  photosphere  and  we  have 


SL 


T 


SL  . 

P 


From  equation  (13),  we  obtain 

ST  SL 

7 —  t(t)  —■ 2  .  (18) 

1  P 

For  a  partial-toroidal  plasma  of  majoi  ladius  R  and  minm  radius  a  (R  a  >> 
1),  the  inductance  is  (Bateman,  1978) 


1 1 


L 


,  4n0R  .  (8R)  „ 

L  =  — rr~  In  —  -  2 
p  2  la  ) 

r  c 


with  =  0. 
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Substituting  the  variation  of  into  equation  (17),  we  find 
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dR  +  ;-l 

(1  dR 

1  da) 
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d© 

R 

dZ  +  L 

R  dZ 

a  dzj 
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dZ 

where  L  =  ln(8R/a)  -  2.  In  order  to  determine  da/dZ,  we  must  related  the 

changes  in  the  pressure  to  changes  in  the  field.  From  equations  (6),  (16) 
and  (17),  we  obtain 


xa  a  cfi  ,-1)1  dR  1  d0  (.  :-lU  da 

66  =  21-0  -  £  1  -*  L  1  +  eL  -Tty 

pi  pj  l  )R  dZ  0  dZ  l  Ja  dZ 


Equating  this  expression  to  the  right  hand  side  of  equation  (9)  and  using 
equations  (10),  (15)  and  (17),  we  can  solve  for  da/dZ.  After  some 
straightforward  algebra,  we  find 


da  _  a_  ‘-1 
dZ  "  2Z  tL 


(l  -  26p)  •  * 


x  2e  1 


-1)Z  dR 


-  z  dR)  Z  d0(„ 

-  YP  R  dzj  +  0  dZ  l2E 


where 


*  -  2 

B  /8rt 
P 

For  the  parameter  values  to  be  used  later,  this  quantity  ranges  from  0.01 
to  0.1.  This  means  that  the  minor  radius  expansion  is  generally  much 
slower  than  the  ;  pex  motion.  Using  these  results  in  equation  (8),  we 
finally  obtain  the  linearized  equation  for  the  height  of  the  apex: 


=  T(t)6Z, 
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where 


1  -  2e  1  - 


««•[’(> 


2c('  -  6p)ldT  ■ 


Here  M  =  rta  nm^  is  the  mass  per  unit  length.  I  f  we  set  z  -  0  (  $__  >>  ♦  ), 
the  right  hand  side  is  positive,  indicating  that  the  perturbation  grows. 
If  we  set  z  =  1  (the  current  is  closed  in  the  photosphere),  the  right  hand 
side  is  negative  so  that  the  displacement  is  restored.  Thus,  there  exists 
a  quantity  e  with 


0  <  z  <1 
cr 

such  that  the  d^(&Z)/dt“'  =  0  for  z  z  . 

el¬ 
and  after  some  algebra,  we  find 


By  setting  F(t=0)  equal  to  zero 


Z  dR  (,  Z  d0  V 
”  ql  R  dZ  ‘  q2  r  0  dZ  nj  J  x 


x  [L-1  -  qj(!  . 


,1Z  dR 

q4 R  dZ  lq 


,)Z  d0 

q4 ^ j  0  dZ  +  q4n 


Here  we  have  defined  h  =  (8rt/B^)  (Z/H)pa ,  q^  =  (1  2(3^)  +  1/2(1  *  2  0^ )  q0 
■  1/2(1  -  20  ),  q3  =  (1  -  20  )“  *  2(1  -  0  )*.  and  q4  s  (1  0  )L  1 .  Note 
that  is  a  function  of  equilibrium  quantities  above  the  photosphere 

only.  A  current  loop  with  z  <  c  is  unstable  to  perturbations  and  a  loop 

.  .  5 

with  z  >  z  is  stable.  For  solar  current  loop  parameters  (e.g.,  Z  -  10 

Sl  -2 

km,  a  ~  10  km,  p  ~  a  few  dynes  cm  '),  z  is  typically  0.2  to  0.4  (Sec. 
III).  The  quantities  e  and  e  have  the  following  physical  interpretation. 
For  z  <  £  <<  1,  the  loop  above  the  pliotosphere  is  a  small  fraction  of  the 

entire  current  distribution.  As  t  lie  loop  expands,  the  changes  in  the  loop 
magnetic  field  and  average  internal  pressure  are  relatively  small  in 
comparison  with  the  changes  in  the  ambient  pressure.  In  paiticulai,  £0^  > 
0  so  that  the  loop  is  unstable.  For  z  >  e  ,  the  loop  is  a  larger  fraction 
of  the  total  current.  The  magnetic  field  and  internal  pressure  decrease 
more  rapidly  in  such  a  way  that  the  displacement  is  restoied.  In  a  more 
formal  sense,  the  present  time-dependent  problem  requires  specification  of 


initial  conditions.  Given  a  loop  above  the  photosphere,  the  quantity 
e( t -0 )  parametrizes  the  initial  configuration  according  to  the  submerged 
current  structure.  This  basic  behavior  depends  on  the  existence  but  not 
the  details  of  the  submerged  current  structure  which  behaves  differently 
from  the  loop  above  (i.e.,  much  less  mobile). 

For  the  unstable  case,  equation  (21)  yields  the  exponential  growth 
time  x  is  given  by 

,-1/2 

x  -  [r(0)J  •  (24) 

It  is  significant  to  note  that  x  «  I  *  so  that  unstable  loops  with  larger 
I  linearly  grows  faster.  For  the  stable  case,  the  loop  can  oscillate 
about  the  equilibrium  position. 

At  this  point,  it  is  useful  to  consider  the  energy  budget  of  a  current 
loop  and  provide  a  more  transparent  meaning  for  the  terms  in  equation  (7). 
Assuming,  for  simplicity,  that  the  loop  is  a  half-torus,  the  total  magnetic 
energy  of  the  semi- toroidal  loop  above  the  photosphere  is  the  sum  of  the 
poloidal  magnetic  energy  and  toroidal  magnetic  energy  where 

E  =  L  I2,  (25) 

p  2  p  t 

where  given  by  equation  (19),  and 


■  8n  ("2A)- 


Using  the  principle  of  virtual  work,  we  find 


rp  -  ^  ['"  IP)  -  >]■ 


k  -  >)• 


where  and  Ff  are  the  major  radial  forces  acting  on  the  entire  loop  due 
to  J^Bp  and  JpB^,  respectively.  It  is  straightforward  to  show  that  the 
total  pressure  force  in  the  major  radial  direction  is 


(29) 


Figure  2  shows  the  various  local  force  components.  Locally,  the  two 
components  of  the  Lorentz  force  are  both  along  the  minor  radius  as  shown. 
However,  when  these  forces  are  integrated  over  the  toroidal  volume,  we  see 
that  J^B  contribution  points  outward  along  the  major  radius  (eq.  [27))  and 
J  pB ^  contribution  points  inward  (eq.  | 28  ] ) .  This  is  entirely  due  to  the 
curvature  of  the  current  distribution.  Adding  the  three  forces  and 
dividing  the  sum  by  rtR  to  get  the  total  force  per  unit  length,  we  recover 
equation  (7),  providing  a  heuristic  derivation.  The  expression  for  F 
shows  that  as  the  major  radius  expands,  the  B  component  does  woik  on  the 
loop,  losing  energy  to  the  loop.  At  the  same  time,  the  loop  does  work  on 
the  Bf  component  so  that  the  Bf  component  gains  energy  as  the  loop  expands. 
Because  the  minor  radius  expands,  the  internal  gas  and  B;  do  work  against 

the  ambient  pressure  and  lose  energy.  On  balance,  there  is  a  net  loss  of 

poloidal  magnetic  energy  to  the  kinetic  energy  of  tiie  loop.  A  fraction  of 
this  energy  is  then  converted  to  thermal  energy  via  drag  heating. 

C.  The  Behavior  of  an  Expanding  Cut  lent  Loop 

In  the  preceding  section,  we  have  described  the  near  equilibrium 
behavior  of  a  model  current  loop  embedded  in  a  background  plasma.  In  this 

section,  we  will  consider  tire  long  time  scaling  behavior  which  will  bo 

useful  for  interpreting  the  numerical  results  to  Ire  obtained. 

As  the  loop  expands,  the  velocity  of  the  apex  increases  and  the  drag 
on  the  ambient  gas  becomes  important.  As  a  simple  model,  we  write 


(viav?)- 


wtiere  F^  is  the  drag  force  pet  unit  length.  V  dZ  dt  is  the  velocity  of 
the  loop  (i.e.,  the  apex),  n  is  the  local  ambient  density  and  <_  ,  is  the 

cl  n 

drag  coefficient.  An  or det  of  magn i t tide  estimate  foi  the  characteristic 
terminal  velocity  in  the  nonlinear  expansion  phase  ran  be  obtained  by 
equating  F^  to  the  driving  force  F^  given  by  equation  (7).  Uc  then  obtain 


V  ~  I  fc,  m.c'n  aR]  ^ 7 ' 
*  t  (  d  l  a  ) 
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1/2 

Equation  (31)  shows  that  V  is  proportional  to  I  /n  •  If  we 

10  9  -3  1  a  5 

estimate  V.  by  taking  I  =  5  x  10  A,  n  =  4  x  10  cm  ,  R  =  10  km,  a  = 

A  ^  a  2  -  1 

10  km  and  using  =  1,  we  find  -  2  x  10  km  sec  .  It  is  of  interest 

to  compare  this  value  to  the  estimated  sound  speed  Cs  in  the  corona: 


C 


s 


1/2 


For  T  ~  2  x  10^  K  and  y  =  5/3,  C  =  2.3  x  10^  km  sec  Although  the 

3  S 

actual  expansion  velocity  depends  on  e,  the  above  comparison  indicates  that 
the  peak  expansion  velocity  can  be  comparable  to  the  sound  speed  under  the 
action  of  curvature  forces  alone.  It  will  turn  out  that  equilibrium  loops 
of  the  type  used  here  can  only  produce  subsonic  expansion.  However,  if  a 
loop  is  allowed  to  be  out  of  equilibrium  initially,  carrying  a  sufficiently 
laLge  (e.g.,  loss  of  equilibrium  at  t  0),  then  it  may  be  driven 

.supersonically  (or  super  Alfvenically  for  magnetized  ambient  plasmas). 
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III.  EVOLUTION  OF  MuuEL  CURRENT  LOOPS 


In  the  preceding  sections,  we  have  discussed  in  detail  the  linear 
(SZ/Z  <<  1)  behavior  of  a  model  current  loop  in  a  background  plasma.  Th" 
description  of  the  long-time  behavior  has  been  limited  to  scaling  laws.  We 
will  now  attempt  to  provide  a  more  quantitative  discussion  of  the  nonlineal 
behavior  by  numerically  integrating  the  equations  of  motion.  Numerical 
examples  are  given  to  illustrate  the  range  of  behavior  under  the  action  of 
curvature  forces  using  parameters  compatible  with  the  solar  environment. 
The  basic  physics,  however,  is  not  limited  to  the  sun. 

As  the  apex  rises  with  increasing  velocity,  equation  (7)  will  no 
longer  be  adequate.  In  addition,  the  drag  due  to  the  ambient  plasma  gas 
may  become  significant.  We  will  incorporate  this  effect  by  a  simple  drag 
model,  equation  (30).  Adding  to  equation  (2),  we  obtain 
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where  M  =  rca‘ p.  Higher  ordei  contributions  are  neglected  foi  simplicity. 
We  have  directly  integrated  the  set  of  equations  (3)  and  (37)  for  a  variety 
of  loop  parameters.  As  stated  in  connection  with  equation  (20),  da/dt  is 
typically  one  tenth  of  dZ/dt  nr  less  so  that  the  minor  radius  is  nearly  in 
equilibrium  for  small  to  moderate  dZ/dt.  We  have  found  equation  (6)  to  be 
nearly  true  even  for  velocity  V  up  to  o.bC^.  This  justifies,  a  poster  Or i , 
the  use  of  equation  (6)  in  the  linear  perturbation  analysis. 

The  drag  term  in  equation  (.32)  is  the  force  which  the  expanding  loop 

experiences  in  displacing  the  ambient  gas.  The  drag  coefficient  is  the 

coupling  coefficient  between  the  loop  and  the  ambient  medium.  In  out 

model,  we  adopt  a  simple  c,  based  on  a  straight  cylinder  transverse  to  the 

flow  in  a  compressible  gas.  Foi  the  subsonic  tegime  with  a  Reynolds  numbei 

R  of  If)6  to  10^,  c,  is  0.5  to  1  (Tritton  1977  ).  At  Mach  1,  c  .  attains  a 
e  d  d 

maximum  value  of  approximately  7  and  decreases  rapidly  foi  latget  Mach 
numbers.  The  supersonic  drag  coefficient  is  obtained  from  Hoernet  (1^31). 

In  our  calculation,  the  ambient  plasma  has  no  background  field.  If 
there  are  ambient  magnetic  fields,  the  drag  coefficient  c ^  must  be  modified 
and,  for  super-Al  f  veil  ic  motion,  MHD  shocks  are  generated.  We  do  not  treat 
shocks  per  se.  The  physical  picture  is  simply  that  if  the  apex  is  driven 
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supersonic  or  super-Al fvenic ,  then  shocks  are  generated.  We  believe  that 
this  treatment  is  a  reasonable  one  unless  the  ambient  fields  are  comparable 
to  or  exceed  the  loop  fields  ( ~20G  for  the  supersonic  examples).  The 
results  are  to  be  interpreted  as  order-of -magni tude  estimates. 

Figure  3(a)  shows  the  velocity  of  the  apex  for  a  loop  with  the  initial 
equilibrium  height  Zq  =  lO^km,  a^  =  2  x  lO^km  and  I  =  4.5  x  10^'  A. 
corresponding  to  =  4 . 5G  and  =  8.1G.  This  loop  is  half  a  torus  (Zq  = 
sq)  with  the  footpoint  separation  of  2so  =  2xlo\m.  This  is  a  case  with 

relatively  weak  magnetic  fields.  The  ambient  pressure  is  taken  to  be  p 

-26  q  ?3 

2  dyn  cm  at  T  =  2  x  10  K  so  that  the  number  density  is  n  =  4  x  10'  cm 

For  this  loop,  we  have  £  =■  0.28  (eq.  [23]).  The  values  of  e 

significantly  smaller  than  e  should  give  rise  to  instability.  Curves  1 

and  2  correspond  to  e(t^O)  -  0.01  and  e(t=0)  =  0.05,  respectively.  The 

2  - 1 

velocity  is  normalized  to  the  sound  speed  =  2.4  x  10  km  sec  .  These 

curves  describe  two  loops  of  apparently  identical  appearance  above  the 

photosphere  with  different  initial  conditions  e(0)  corresponding  to 

different  submerged  structures.  For  Curve  1,  the  flux  enclosed  by  the 

entile  current  distribution  is  one  hundred  times  what  is  above  the 

photosphere  and  for  Curve  2,  the  total  flux  is  20  times  what  is  above. 

Because  of  the  low  current  and  weak  magnetic  field,  these  loops  do  not 

expand  rapidly.  Although  not  shown  here,  these  loops  continue  to  expand 

slowly  even  after  one  hour  with  the  major  radius  reaching  1.5  to  2  times 

the  initial  values.  The  expansion  is  nearly  exponential  for  the  first  20 

minutes.  In  Figure  3(b),  t he  height  of  the  loop  apex  is  show.  In  general, 

with  other  parameters  being  equal,  loops  with  smaller  values  of  t  <  £ 

expand  more  rapidly  expand  more  rapidly  to  larger  values  of  Z,  and  in  cases 

where  loops  can  attain  "second"  equilibrium  (see  Fig.  4  for  an  example), 

they  do  so  later  and  at  larger  values  of  Z.  Also,  as  a  loop  expands,  the 

expansion  tends  to  slow  down  because  the  current  and  magnetic  field 

decrease  and  e ( t )  increases,  sometimes  reaching  a  second  equilibrium.  The 

dashed  line  describes  a(t)/'a  ,  showing  that  the  minor  radius  expansion  is 

much  less  than  the  increase  in  the  apex  height. 

Figures  4  shows  the  behavior  of  a  smaller  loop  with  an  initial  height 

of  Zn  -  lO^km  and  R  -  7.2xl0^km  with  Z  /s  -  1.5.  The  footpoint 

U  °  4  on  ' 

separation  is  2.s  =  1.33x10  km  with  0/n  -  0.374.  The  aspect  ratio  is  taken 

°  3  q 

to  be  5  so  that  an  =  1.4  x  10  km.  The  current  is  I  =  3.3  x  10  A  with  B 
0  t  p 

4.5G  and  B^  =  8.1G.  For  this  loop,  we  find  £  =  0.05.  In  this  example, 

we  have  used  £  =  0.03  and  the  loop  is  only  mildly  unstable.  The  apex 
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velocity  reaches  a  maximum  of  ~0.  1C  with  a  rise  time  of  7  minutes.  The 

5  4 

velocity  then  decreases  until  the  height  reaches  Z  =  6.5x10  km. 

Subsequently,  the  apex  executes  damped  oscillation  with  a  period  of  -5  min. 

Preceding  examples  show  model  loops  which  are  fairly  "tali",  i.e., 

with  Z  >  s  ;  Z  =  s  for  Figure  1  and  Z  -  1.5s  for  Figure  4.  Flatter 

equilibrium  loops  (Z  <  s  )  are  also  easy  to  find.  It  seems  intuitively 

reasonable  (perhaps  necessary)  that  flatter  loops  should  be  able  to  slowly 

evolve  to  taller  configurations  if  small  changes  in  parameters  aie  made. 

This  is  in  fact  the  case.  In  Figure  5,  we  use  a  loop  with  Z  -  lu  km,  s 
4  4  no 

1.25x10  km,  R  =  1.28x10  km,  R  /a  =5  and  £(t=0)  =  0.05.  The  cur  lent  :s 
g  0  o  o 

I.  =  6x10  A  with  B  -  4.50  and  B  =  8.10  in  equilibrium.  At  t.n,  We 

t  p  t  1 

increase  the  current  I  ftom  the  above  equilibrium  value  by  1 % .  This,  of 

course,  corresponds  to  increasing  the  magnetic  field  twist  slightly.  (Note 

that  the  eqs.  (3)  and  (32)  ate  not  limited  to  equilibrium.)  The  result 

shows  that  the  apex  rises  and  executes  damned  oscillation  about 

4 

1.16x10  km.  The  expansion  velocity  is  small,  nevet  exceeding  ~0. 10  . 
Analogously,  if  the  current  is  decreased  slightly,  the  height  can  dectoase 
slowly  followed  by  small  amplitude  damped  oscillation.  This  behavior 
suggests  that  a  current  loop  can  evolve  to  a  taller  (lower  )  loop  in  a 
quasi  equilibrium  manner  if  the  current  increases  (decreased)  slowly. 

In  general,  with  all  other  qtian  t  i  t  i  ec:  being  equal,  smaller  unstable 
loops,  have  shorter  e  folding  times  (eq.  |21|)  because  of  the  reduced 
inertia.  Unstable  loops  with  latget  currents  I  al-o  have  shorter  e 
folding  times  because  ot  the  increased  Ioretur  tot  re.  Iti  figure  t. ,  we  show 
an  example  with  latget  cur  rents.  The  loop  is  not  in  equilibrium  initially. 
This  example  may  be  relevant  to  a  loop  which  suffers  loss  of  equilibrium. 
The  parameters  used  are  Z,.  -  in  km.  7  -  s  .  a,  -  lO^km  and  1  10^A  so 

that  B  -  20G  and  B{  -  210.  Curve  1  corresponds  to  £(0)  -  tt.Ol  and  the 

loop  attains  Mach  3  in  less  titan  one  minute.  At  t  =  18  min,  the  apex 
velocity  is  approximately  V/C  -  5  or  V  1200  km  s  1  and  the  apex  has 

risen  to  Z  -  20  x  In’3  km.  Subsequent!',',  the  expansion  volo<  ity  slowly 
decreases  over  tens  of  minutes  as  the  loop  expands.  Curve  2  mi  responds  to 
e(0)  =  0.05  and  the  configuration  is  slower  than  l ha l  desci i bed  by  Cm  ve  1. 
The  apex  attains  Mach  3,  the  maximum  vr  lor  ity,  in  about  5  minute'  .  also 
with  a  sharp  rise  neat  t  o.  The  velocity  « lien  slowly  dec  i  eases  in  Mach 
in  about  10  minutes  as  the  loop  expands.  For  smaller  cut  rents,  the 
velocities  are  smaller.  Ftom  Figures  1  -  6,  it  is  clear  that  curvature 

force  can  produce  a  wide  range  of  behavior. 


IV.  PHYSICAL  IMPLICATIONS  AND  DISCUSSION 


We  have  described  the  dynamics  of  the  apex  of  a  model  current  loop 
embedded  in  a  background  plasma.  The  structure  is  such  that  the  semi- 
toroidal  section  of  the  loop  is  in  the  upper  tenuous  plasma  while  the 
remainder  of  the  current  distribution  is  embedded  in  a  much  denser  plasma. 
The  dynamical  properties  obtained  are  most  applicable  to  the  apex  of  the 
semi- toroidal  loop.  We  reiterate  that  the  objective  is  not  to  model 
specific  phenomena  but  to  understand  the  basic  effects  of  curvature  forces. 
For  this  purpose,  we  have  constructed  the  model  in  such  a  way  that  the 
model  loop  behavior  is  primarily  determined  by  the  curvature  forces.  In 
the  preceding  sections,  we  have  given  theoretical  and  numerical  results. 
Although  the  results  cannot  be  applied  directly  to  observed  solar  phenomena 
because  of  simplifying  assumptions,  it  is  useful  to  determine  the  range  of 
expected  behavior  using  parameters  compatible  with  the  solar  environment. 

ohservat ionally,  it  is  not  always  easy  to  determine  the  magnetic 
stinitnie  or  its  motion.  However,  signatures  of  motion  may  he  manifested 
as  heating  of  coronal  gas  and  moving  gaseous  material.  Here,  we  will 
examine  some  possible  observational  implications.  For  this  purpose,  it  is 
useful  to  consider  the  rate  at  which  the  magnetic  energy  is  converted  to 
thermal  energy  via  drag.  Wc  have  calculated  the  quantity 


fot  the  model  loops  described  in  the  preceding  section.  Heie,  F^  is  the 
drag  force  given  by  equation  (30)  and  dE/dt  is  the  rate  at  which  the 
ambient  gas  is  heated  by  drag  due  to  the  apex  motion.  In  calculating  this 
quantity,  we  have  assumed  that  only  one  third  of  the  semi -torus  around  the 
apex  is  effective  in  drag  heating.  As  the  above  expression  indicates,  the 
heating  rate  is  proportional  to  VJ.  We  have  also  computed  the  time 
integrated  total  energy  which  the  magnetic  field  loses  in  the  form  of  loop 
plasma  kinetic  energy  and  drag  heating.  This  quantity  is  essentially  equal 
to  the  time-integral  of  dE/dt  plus  the  loop  kinetic  energy.  As  pointed  out 
before,  the  minor  radial  expansion  is  found  to  be  much  slower  than  the  apex 
motion  so  that  it  is  neglected  in  comparison  with  the  majoi  radial 
expans i on . 

Figure  7  shows  tfie  energy  release  rates  due  to  drag  for  the  loop 

described  in  Figure  3.  For  £  -  0.01  (Curve  1),  the  rate  reaches  lO^  erg 

1  2  G  1 

s  at  t  =  20  min  and  increases  to  10  erg  s  at  t  =  30  min.  During  this 
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time,  the  major  radius  increases  from  1  .  1  5Z()  to  l.oRy  Before  t  20  min, 

the  loop  exhibits  only  slow  motion  and  insignificant  energy  output.  For  £ 

=  0.05  (Curve  2),  the  loop  motion  is  less  pronounced  with  the  energy 

25  l 

release  rate  in  the  range  of  10*'  erg  s  during  t  -  20  min  to  t  =  10  min. 

Figure  7(b)  shows  the  t ime- integi ated  energy  converted  from  the  magnetic 

field  to  thermal  and  kinetic  energv.  For  e  -  0.01  (Curve  1),  the  total 

29 

amount  of  magnetic  energy  teleased  in  .10  minutes  's  ~  1.5  x  10  "  erg  while, 

for  e  =  0.05  (Cutve  2),  it  is  10*  erg.  For  both  cases,  roughly  one  half 

of  the  energy  is  in  the  form  of  thermal  energy.  In  Figure  R.  we  give  the 

energy  output  protile  tot  the  loop  desci i bed  in  Figute  6.  For  this  loop, 

the  magnetic  tield  components  ate  -  .'Ot,  and  the  apex  can  he  driven 

supersonic  with  correspondingly  greater  magnetic  energy  i e lease .  Howevet  , 

as  discussed  before,  this  loop  is  not  initially  in  equilibrium  in  t  lie- 

context  of  the  present  model.  'For  curve  1  (s  -  o.oi),  the  maximum  eneigy 

release  rate  is  roughly  •*  In*  ^  t  i  g  *  with  a  time  integi  ated  total  ol  '• 
1 1 

x  In  erg  in  It '  minute;  .  Curve  7  (f  o.l).  shows  an  energy  output 
profile  in  which  t  lit  peak  heating  <>.  <  in  s  in  ,i  dutat  ion  ot  l11  minutes  with  a 

long  decay  phase  lusting  for  ten*  ot  minutes.  Th*  total  emu  gv  released  is 

>,  i 

roughly  3  x  H>  erg  in  »>  minutes.  F<o  » 1 1 <  •; c-  r  ui  ves,  their  is  a 

possibility  of  strong  shock  heating.  No'e  • bar  magnet  if  energy  (onveisjon 
takes  place  on  t  he  * ime  s<ule  o(  tens  of  -n  mites.  The  rate  of  energy 
release  cai>  have  a  wide  range.  This  type  of  "dynami  <-a  1 "  magnet  i<  eneigy 
release  mechanism  has  hern  suggested  ear  1  ier  (Xue  ami  Chen .  r>8<>) .  As 
mentioned  eat  1  iei  ,  the  impoitame  o  1  mec  ban  i  <  a  1  energy  output  in  t  he  flair 
energy  budget  has  been  discussed  ( Veldt  et  a  1  .  .  l'*HO) .  Also,  it  a  1  a  t  g< 
number  of  loops  release  eneigy  at  slow  tales,  they  may  <  ont  r  i  tune  to  non 
violent  theimal  energy  input  to  the  corona. 

It  is  of  interest  to  estimate  the  temperature  of  the  ambient  gas  vh i '  1. 
is  heated  by  the  supersonic  motion  of  the  apex.  Foi  strong-  shocks,  the 
temperature  T+  behind  the  slunk  front  can  he  estimated  by 


(handau  and  Lifshitz,  1939)  where  M  i  the  Hath  number  of  i  he  shod  and  I 
is  the  ambient  t  empera  t  ui  e .  Taking  M  -  1  (Fig.  6(a))  and  y  ..  3/3,  we  find 

T^  =  3.7  T  .  Using  2  x  10'’  K,  we  find  7.4  x  1  o*’  K.  For  larger 

values  of  M,  the  temperature  is  higher  •  Thus,  in  this  particular  example 


(Curve  2),  the  coronal  gas  in  the  vicinity  of  the  apex  could  be  heated  to 
approximately  10  K  and  the  heated  gas  might  be  seen  to  travel  away  from 
the  sun  with  a  peak  value  of  ~  800  km  sec  This  phase  can  last  for  tens 
of  minutes  with  the  velocity  and  heating  diminishing  with  time.  For  Curve 
1,  the  velocity  is  considerably  higher  (a  peak  value  of  ~  1200  km  s  ^). 
These  calculated  apex  velocities  are  similar  to  those  associated  with 
certain  dynamical  effects.  For  example,  the  sources  of  moving  type  IV 
bursts  are  quoted  to  have  velocities  in  the  range  of  roughly  200  km  s  ^  to 
1500  km  s  ^  with  typical  velocities  around  400  km  s  ^  (e.g.,  Svestka, 
1980).  For  coronal  mass  ejection  events,  Gosling  et  al.  (1976),  fot 
example,  reported  speeds  ranging  from  less  than  100  km  s  ^  to  more  than 
1200  km  s  *  with  the  average  being  roughly  500  km  s  ^ .  They  also  reported 
association  of  type  II  radio  bursts  with  events  moving  faster  than  about 
400  km  s  Note  that  we  do  nor  attempt  to  "explain"  these  phenomena  here. 

Mouschovias  and  Poland  (1978)  and  Anzet  (1078)  have  described  loop- 
type  (as  opposed  to  "bubble")  fiansient  models.  in  anzet's  work,  poloidal 
cm  rent  density  J  and  plasma  pressure  are  negiected,  effectively  treating 
the  loop  as  a  metallic  viie.  The  loops  are  not  not  in  equilibrium 
initially.  Ir.  addition,  the  ambient  plasma  is  neglected.  Nevertheless, 
t lie  underlying  physics  is  similar  to  that  of  out  model  in  that  both  models 
use  Loientz  curvature  forces  to  drive  current  loops.  Hete.  although  the 
present  model  is  not  specifically  for  coronal  transients,  we  make  a 
qual i tat i ve  comparison  with  the  work  of  Anzet .  In  Anzet's  wot  k .  it  was 
found  that  magnetic  fields  of  1G  can  drive  coronal  transients.  In  out 
model,  ve  estimate  the  necessary  magnetic  fields  to  tie  greater.  This  can 
he  unde! stood  in  the  following  way.  In  Anzet's  model,  the  poloidal  current 
density  .1  and  pressure  gradient  are  neglected.  In  the  toroidal  geometry, 
the  fot<e  J  B(  acts  to  counter  the  expansion  of  the  apex.  Furthermore,  the 
ambient  material  which  also  acts  to  oppose  ihe  expansion  is  neglected. 
The  only  retarding  force  is  gravity.  In  our  present  model,  ve  include  the 
poloidal  current  and  plasma  pi essiii e.  Fot  the  nonequ i 1 i hi i urn  examples 
(Figs.  6  and  8),  the  magnetic  !  i'  Id  components  ate  It1  .’(>G  in  the  level 
cot ona  so  that  gravity  is  unimportant  (see  below).  Thus.  Anzet's  model 
tend  to  requite  smaller  magnetic  fields  than  out  model  to  diive  cut  tent, 
loops  to  a  given  velocity.  In  addition.  the  (intent  loops  used  bv  Anzet 
ate  much  larger,  initially  -  n.5R  .  In  nut  model,  the  magnetic  field  is 
also  weaker  at  comparable  altitutes.  Taking,  for  example,  Curve  1  of 
Figttie  6,  ve  find  that  at  t  10  min,  Z  -  IR  ,  the  magnetic  field 


c omponen  t  s  ate  i  oughly  lot;.  Fm  Cutve  ?  at  t  10  min,  Z  =  ,’R^ .  the 

magnetic  tield  components  ate  roughly  fit;.  We  attribute  the  numerical 
d!,retences  piimatily  to  the  neglect  ot  poloidal  current  and  plasma 
pressure.  In  addition,  the  inclusion  of  ambient  coronal  gas  allows 
conversion  of  magnetic  energy  to  thermal  energy  in  our  model. 

Citing  the  tesult  of  Anzer  t  1978).  Carlqvist  and  Alfven  (1980) 
suggested  that  expanding  current  loops  may  contribute  to  the  acceleration 
of  solar  wind  particles.  By  viitue  of  including  both  toroidal  and  poloidal 
current  components,  plasma  ptessuie  and  ambient  gas,  out  model  can  be 
extended  to  give  a  potentially  better  description  of  loop  s true  true  and 
dynamics  beyond  the  immediate  vicinity  of  the  sun.  (Also  include  eq •  [  hi] 

if  gravity  is  important.)  We  note  that  Klein  and  But  laga  (198?)  reported 
observational  evidence  ot  magnetic  clouds  at  1  All  consistent  with  the 
geometry  ot  magnetic  '.oops.  f  t  lias  been  suggested  (t'occoni  et  al .  .  1958; 
Cold  1959,  19A?)  that  the  magnetic  field  Litres  of  magnetic  clouds  may  bo 
anchored  in  the  sun.  However,  t  lie  <‘c'  l!-’1  magnetic  'tincture  and  dynami<  s 

of  magnetic  clouds  in  the  m t f t p ■  niotat  •  pair  has  not  been  lullv 
understood.  It  is  possible  thit  expa  >  c  uiten  f  < .  i  magnet  i<  loops  can 

na  t  urn  1  I  i  •  e  t  i  sr  t  o  t  lice.  tine  i  it  '  ■  .  it  i  si  gu  if  i  ■  ant  to  note  that 

Klein  anti  Bur  la  J  i  (19B?)  suggested  1  ha  t  cei  tain  diffeient  types  of  •  louds 
•nay  be  different  man '  m-  ;  >  a  t  i  uin.  t  inglt  phenomenon  such  as  coronal 

t  [  aiis  t  on  t  s  . 

The  mot  ion  |o  !  |r.  ;  ■  i i ■ .  ■  i-  t  in  ha  '■  o|  t  !,<■  ■  nnvcct  ion  cone  ha'  been 

mr  I'Ihi  I  ( t'hn'i  md  1  ;  h.  :  ,  '  ;  !' :  In  :  rn-i  1  'in  hi  .  19,8 7  )  .  it  mi.  Ii  1  1  u:- 

tubes  ,u  c  am  Inn  nd  bein'.  '  in  un  •  t  ion  con.',  t  hen  iqu  i  1  i  hi  i  uv  and 

dynam  i  r  a  I  <  mt-  id.ua>  ion-  i-  -  i  i  h.  :  in  a  f  and  in  cue  and  flten  (  1  9  get;  ]  »H  •  ) 

r.av  I't  app  1  i  i.  a  b  1 1  .  In  i  • .  1 1  t  i .  a  t  i .  i :  .  t  In  i  nd  m  t  i  •.  i  « •  i  t  (•<.  t  ■:  des.  i  i  hod  by  f  may 

% 

he  i  (leant  . 

In  tin-  examples  .  i .  t  hi  ■  pap.  i  ,  tire  i  ole  uf  gtavity  has  not 

Imen  t  ons  i  dot  cd  .  Fm  ome  j  ?  1 1  •  r  i .  >  :v .  no  .  gi  .r  it.  ;->,iy  be  i '"peu  t  an  *  .  lot  t  lie 

apex  ot  a  loop.  gta  ii,  ,n  i  a  lung  tin  mm  i  .  i  <  1  i  n-  , .  t  la  •  i 

s  t  i  a  i  gh  t  f  ot  va  id  to  >  in  I  nd<  tin  gi.rit.itiuii.il  t  ni  .  I  .  In  a  . 
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If  n  >  n  (e.g. ,  coronal  mass  ejections),  F„  is  downward.  li  n  <  n  ,  then 
a.  u  a 

the  structure  is  buoyant  and  F^  is  upward.  For  magnetic  fields  of  10 

20G,  the  curvature  forces  dominate  the  gravitational  force.  ,  For  example, 

-9  -3 

for  the  supersonic  loop  depicted  in  Figure  6  wi th  a  density  of  10  cm  ,  the 

curvature  forces  are  of  the  order  of  lO^dyn  cm  ^  while  F_  is  of  the  ordet 

8  1  b  /  *•' 
of  10  dyn  cm  .  The  basic  tenets  of  curvature  effects  remain  qualitatively 

valid  with  the  addition  of  gravity. 

In  summary,  we  have  theoretically  studied  the  behavior  of  a  simple 
semi  toroidal  current  loop  under  the  action  of  curvature  forces.  It  has 
been  shown  that  such  loops  are  capable  of  exhibiting  a  wide  range  of 
dynamical  behavior.  Starting  with  MHD  equilibria  with  the  curvature  forces 
explicitly  balanced  and  with  the  footpoints  of  the  loop  remaining 
stationary,  Loops  can  expand  with  a  wide  range  of  subsonic  velocities, 
giving  rise  to  a  correspondingly  wide  range  of  magnetic  energy  output. 
Some  loops  can  attain  second  equilibria.  Loops  may  also  evolve,  in  a 

quasi  equilibrium  manner,  from  "flat"  to  "tall"  loops  if  the  current  is 
increased  slowly.  The  typical  time  scales  for  motion  and  energy  release 
are  tens  of  minutes.  Given  loops  initially  in  equilibrium  with  no  ambient 
magnet i'-  fields,  the  subsequent  motion  seems  to  be  subsonic  with  relatively 
slow  heating  of  the  coronal  gas.  If  ve  start  with  nonequi 1 ibtium  loops 
with  laige  currents,  possibly  as  a  result  of  loss  of  equilibrium,  they  can 
attain  highly  supersonic  (or  supet  Alfvenic  in  magnetized  ambient  plasmas) 
ex  pa  i  is i on  velocities  with  rapid  heating  due  to  shock  heating.  Thus,  it 
appeals  that  expansion  of  magnetic  loops  can  contribute  energy  to  the 
corona  I  gas.  along  with  other  possible  mechanisms,  with  a  wide  range  of 
enetgv  release  rates. 

A  novel  but  somewhat  unconventional  featute  of  the  model  is  the 
inclusion  of  submerged  current  distributions  in  the  dynamics  of  the  loop. 
As  pointed  out  previously,  the  present  model  depends  on  the  existence  but 
not  any  details  of  tire  submerged  cut  rent  distributions.  The  submerged 
mi t out  is  parametrized  by  the  quantity  £,  the  latio  of  fluxes,  defined  by 
equation  (14).  Although  not  measurable  in  reality,  this  is  a  physically 
meaningful  quantity.  Tire  condition  for  instability  is  found  to  he  given  hv 
£  <  £'  |  .  It  is  of  interest  to  speculate  on  how  such  an  equilibt  ium  loop 
can  lx  manufactured  by  the  sun  since  a  loop  which  is  observed  to  exist  f  oi 
an  extended  period  of  time  (perhaps  a  few  days)  is  presumably  in  stable 
equilibiium.  For  example,  if  a  flux  loop  merges  through  fhe  photosphere, 
it  expands  until  it  reaches  an  equilibrium  for  which  s  >  e  .  Recall  that 


14 


e  depends  only  on  physical  parameters  above  the  photosphere  while  £  - 
4>^/ .  It  is  possible  within  the  context  of  this  model  that  a  stable  loop 
can  become  unstable  if  increases  (smal'er  e).  This  may  occur  due  to 

topological  changes  in  t  Ire  subphotospherie  magnetic  flux  structure  to  which 
the  loop  is  connected.  There  need  not  be  significant  changes  in  the  loop 
parameters  or  coronal  conditions  prior  to  the  onset  of  expansion. 

In  out  model,  the  current  loop  above  the  photosphere  is  connected  to 
the  submerged  structure  via  flux  tubes  going  through  the  photosphere.  The 
flux  tubes  serve  as  a  conduit  foi  electromagnetic  and  other  processes.  Ve 
have  not  addressed  issues  concerning  the  details  of  possible  flux 
structures  below  the  phtosphere  and  possible  dynamics  of  and  transport 
mechanisms  in  such  current  or  flux  structures  in  subphotospherie  regions. 
It  appears  that  these  properties  can  influence  the  dynamics  of  the  current 
loop  above.  Thus,  in  a  mote  complete  model,  the  parameter’  e  will  contain 
much  mote  physical  effects.  In  our  calculation,  we  kept  only  one  aspect, 
"iz.,  the  relative  t lux  and  the  inductive  effects,  to  i 1  lust  1  ate  certain 
basic  physical  properties.  An  adequate  discussion  of  these  issues  requires 
more  understanding  ot  subphotospher  i c.  fields  and  plasma  properties,  both 
inside  and  outside  of  possible  flux  structures.  As  a  first  step,  ve  have 
included  the  submerged  current  structure  only  as  a  parameter. 
Nevertheless,  the  basic  conclusion  that  the  behavior  of  current  loops  can 
depend  on  submeigeu  ament  structures  is  not  based  on  detailed  assumptions 
and  seems  worthy  of  more  detailed  consider  at  ion. 

As  a  starting  point  for  studying  the  effect;  of  curvature  forces,  ve 
have  used  a  simple  mode  1  system  in  the  limit  ot  veto  background  magnet  i * 
field.  Although  these  equilibria  (  (S  <  <*)  may  be  relevant  to  certain 

structures  sin  h  as  "tool  loops"  with  a  piesmite  deficit  inside  (e.g., 
Foukal  197b;  (.hinder  i  et  al.,  1 '  *  7  7;  Mood  and  Priest  many  of  'he 

observed  loops  are  likely  to  cortespond  to  the  0  >  ('  case.  In  an  upcoming 

paper,  we  will  report  on  loop  equilibria  with  nonwetn  hackgtonnd  magnet  it 
fields  such  that  6  >  ".  The  present  analv-  i  •;  < an  be  applied  to  chi;  :  as* 
in  a  similar,  way.  one  di  f  lot  enee  is  that  gioater  amounts  of  current  can  he 
suppoited  in  equilibrium  and  a  wider  range  of  expansion  morion  can  lie 
exhibited  by  equilibrium  loops.  The  basic  effect;  of  cinvarme  hm  c 
not  limited  to  the  particular  class  ol  equilibria  used  in  thi1*  paper. 
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Throughout  the  paper,  the  emphasis  and  the  objective  have  been  to 
isolate  the  essential  physics  of  curvature  forces  and  to  describe  the  range 
of  possible  behavior  that  may  be  exhibited  by  current  loops.  The  results 
are  to  be  taken  in  an  order-of-magni tude  sense.  The  effects  of  various 
approximations  and  different  types  of  equilibria  should  be  studied  for  a 
more  complete  understanding.  In  this  paper,  ve  do  not  claim  to  explain  any 
particular  observed  phenomena.  The  precise  role  and  the  relative 
importance  or  unimportance  of  curvature  forces  in  any  observed  phenomenon 
must  be  evaluated  in  future  research.  Our  reference  to  observations  is 
limited  to  pointing  out  the  possible  relevance  to  the  sorar  environment. 
Nevertheless,  it  appears  that  current  loops  under  the  action  of  curvature 
forces  can  exhibit  a  range  of  motion  (velocities,  energy  dissipation,  etc.) 
in  the  right  "ball  park"  in  the  '■'cr.iext  of  dynamical  solar  phenomena. 
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CORONA 


Schematic  drawing  of  a  model  current  loop  embedded  in  the  corona. 
Components  of  the  current  density  J  and  magnetic  field  are  shown. 
The  subscripts  "t"  and  "p"  refer  to  the  toroidal  and  poloidai 
directions,  respectively.  The  radius  of  curvature  is  R  and  the 
apex  height  from  the  photosphere  is  Z.  The  footpoint  separation 
is  2sq.  The  segment  between  the  dashed  lines  is  the  "apex 
region",  taken  to  be  1/3  of  the  loop  for  numerical  calculations 
(Secs.  Ill  and  IV).  If  the  apex  region  were  devided  into  thin 
slices,  the  center  of  mass  of  each  slice  would  lie  on  the  dash  dot 
line.  No  particular  structure  is  specified  below  the  photosphere. 
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Behavior  of  a  model  loop  initially  in  equilibrium  with  Zq  =  10  km, 
Z  =  s  and  a  =  2  x  104km.  e  =  0.28  (eq.  (231).  I  = 

O  O,  Cl  1  1  t 

4.5x10  =  4.5G  and  ==  8.1G.  For  both  figures,  Curve  1  is 

e(t=0)  =  0.01  and  Curve  2  is  e(t=0)  =  0.05.  (a)  Velocity  profile 

2  1 

normalized  to  the  sound  speed  C^  =  2.4  x  10  km  sec  .  (b)  Apex 

height.  The  dashed  line  shows  a(t)/a  corresponding  to  Curve  1. 


2  (104  Km)  V/C 


TIME  (MINUTES) 

4 

Fig.  4  Behavior  of  a  model  loop  initially  in  equilibrium  with  Aq  =  10  km, 

s  =  Z  / 1 . 5  and  a  =  1.4xl03km.  I  =  3.3xl09A,  B  =  4.5G  and  -- 
ooo  t  p  t 

8.1G.  =  0.05  and  e(t=0)  =  0.03.  (a)  Velocity  profile.  (b) 

Apex  height. 
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Height  of  the  apex  of  a  quasi -equi  1  i bi  i um  model  loop  vitli  Z  -  l11' 
km,  Z  =  0.8s  with  s  --  I  .  ZbxlO^km  and  a  -  2.6xlo\-m.  I 

gO  O  O  O  t 

6xl0'A,  B  =  4 . 5G  and  B  -  8.10.  The  cm  ve  corresponds  to  e(t-<M 


dE/dt  (1025  erg  sec 


Fig. 
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Magnetic  energy  released  by 
and  a  =  2  x  10*  km).  Curve 
(a)  Rate  of  drag  heating 
magnetic  energy  released  as 
heating  is  roughly  one  half 


the  model  loop  of  Fig.  1  (k  10J  km 
1  is  c.  -  0.01  and  Curve  ?  is  c  0.05. 
near  the  apex.  (bN  Time  integrated 
drag  heating  and  kinetic  energy.  Otag 
of  the  total  energy  released. 
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TIME  (MINUTES) 

Magnetic  energy  leleased  by  the  nonequilibrium  loop  of  Fig.  6  (R  - 
10  km  and  a  =  2  x  10^  km).  Curve  1  is  t  =  0.01  and  Curve  ?  is  z 

-  0.01.  (a)  Rate  of  drag  heating  near  the  apex.  (b)  Time 
integrated  magnetic  energy  released  as  drag  heating  and  kinetic 


energy . 


